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^ ■ ABSTRACT 
O 

We investigate 35 pre-stellar cores and 36 proto-stellar cores in the Perseus 
! molecular cloud. We find a very tight correlation between the physical parameters 

C/^ ! describing the N2H"'' and NH3 gas. Both the velocity centroids and the line 

^ ! widths of N2H"*" and NH3 correlate much better than either species correlates 

with CO, as expected if the nitrogen-bearing species are probing primarily the 
dense core gas where the CO has been depleted. We also find a tight correlation 
'ES ! in the inferred abundance ratio between N2H"'" and para-NHa across all cores, 

with A^(p-NH3)/iV(N2H+)= 22 ± 10. We find a mild correlation between NH3 
(and N2H"*") column density and the (sub) millimeter dust continuum derived H2 
^ [ column density for pre-stellar cores, iV(p-NH3)/A^(H2) ~ 10~^, but do noi find a 

00 I fixed ratio for proto-stellar cores. 

^ I The observations suggest that in the Perseus molecular cloud the formation 

and destruction mechanisms for the two nitrogen-bearing species are similar, 
O I regardless of the physical conditions in the dense core gas. While the equiva- 

lence of N2H+ and NH3 as powerful tracers of dense gas is validated, the lack of 
correspondence between these species and the (sub)millimeter dust continuum 
I observations for proto-stellar cores is disconcerting and presently unexplained. 
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Introduction 



The study of dense cores in nearby clouds like Taurus and Perseus has provided a main 
observational input for testing theories of star formation. These cores are currently forming 
solar-type stars, and the study of their internal structure offers a unique opportunity to 
explore observationally the initial conditions of stellar birth and the first stages of proto- 
stellar evolution. 

The systematic improvement in sensitivity and angular resolution of IR and radio ob- 
servations has gradually revealed that dense cores have a more complex internal structure 
than initially thought. Studies of the density in pre-stellar cores, for example, have found 
radial profiles with a central flattening reminiscent of the isothermal (Bonnor-Ebert) mod- 
els, although the overall state of equilibrium of these c ores is still a matter of debate (e.g.. 



Di Francesco et al 1 120071 : IWard-Thompson et al 1 120071 ). The chemical composition of the 



core material also presents systematic variations with radius. The innermost gas in a core 
prior to star formation is usually depleted of C-bearing molecules like CO, CS, or HCO"'", 



(Caselli et al 


1999; 


Bergin et al 


2002; 


Tafalla et al 


2002) 



cm 

affects our ability to determine the physical conditions and kinematics of the pre-stellar gas, 
as it limits the number of tracers available to observers and adds a layer of complexity when 
comparing the emission from different molecular species. Because of this, understanding and 
characterizing mo lecular freeze out has b ecome a necessary step to realize the potential of 
core studies (e.g. 



Bergin fc Tafallall2007f ) 



Most of our information on molecular freeze out comes from the detailed study of a 
selected group of cores and globules, like L1544, B68, L1498, and L1517B (see previous 
references). These works have shown that under freeze out conditions, the most reliable 
tracers of the dense core material are nitrogen-bearing molecules like N2H"'" and NH3 together 
with the dust component. These three tracers systematically provide consistent views of the 
dense cores, in terms of their maps having similar shape, size, and peak position, and their 
molecular spectra presenting similar central velocity and linewidth. Such a good agreement 
between tracers suggests that their emission arises from the same gas in a core, and that 
their chemical similarities dominate over a number of significant differences between their 
emission properties. T he J, i^=(l, l) transition of NH3, for example, has a rather low critical 
density of ~ 10^ cm~^ fjEvanslll999), while the commonly- u sed J=l-0 transition of N2H"'" has 



a much higher value (~ 10^ cm ^, lUngerechts et al 1119971 ). and the dust emission is directly 



proportional to th e gas column density but sensitive to the dust temperature and emissivity 
faildebrandlll983l ). 



Although the main observational characteristics of molecular freeze out seem well estab- 
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lished by now, a number of unsolved issues require further investigation. Most previous work 
for example, used reduced and highly se lected samples of targets, so detailed radiative trans- 



fer rn odeling could be carried out (e.g., iCaselli et al Ill999l : iBergin et al Il2002t iTafalla et al 
20041 ). Such a strong target selection, unfortunately limits the statistical significance of 



the work and precludes the investigation of topics like cloud-wide chemical variations, time 
evolution of the core composition, and a critical inter-comparison between the behavior of 
the freeze out resistant tracers (N2H"'", NH3, dust continuum). To investigate these issues, 
it is necessary to carry out systematic, multitracer observations of a large sample of cores 
in a cloud, something that was not possible just a few years ago. Fortunately, a number 
of surveys of nearby molecular clouds have been recent ly undertaken acr oss a wide range 
of wavelengths and molecular species (see for example, iRidge et al II2006I ). and these have 
provided unique data sets of continuum emission and molecular line diagnostics. 

One star- formin g region which has received significant recent attention is the Perseus 



molecular cloud (e.g., Hatchell et al 112005 Enoch et al 1 20061: J0rgensen et al 



1 



2006a 



Kirk. Johnstone, fc Di 



2006:lKirk. Johnstone. &: Tafalla 



200 



RebuU et al 



2007 



Rosolowsky et al 



20081: 



Hatchell fc Dunham 



2OO9I ) . This cloud contains clustered low and intermediate mass proto-stellar candidates, and 
seems to represent a case in between the low-mass star-forming cloud of Taurus and the more 
massive Orion star-formi ng region. A dista nce of 250 ± 50 pc to Perseu s has be e n ado pted 



by the Spitzer c2d t ea m (lEvans et al 



20031). follow ing measurements by ICernis I (jl993[ ) and 



Belikov et al.l (120021 ). Ijorgensen et al I ( l2007l . |2008|) combined the Spitzer observations and 
the submillimeter continuum data of this cloud to produce a complete sample of deeply 
embedded protostars in Perseus and to determine the clustering properties of the pre-stellar 
cores an d protostars. From a comb i ned N 2H^ and C^^O survey of this Perseus core pop- 
ulation, iKirk. Johnstone, fc Tafalla I (120071 ) showed that the gas motions in the vicinity of 
submillimeter cores vary significantly, from relatively quiescent inside the individual cores to 
dynamic in the surroundin g gas. Further ob s ervati ons of the Perseus co re population have 



been recently prese nted by iRosolowsky et al I (120081 ). These authors (also Foster et al 112009 



Schnee et al 1120091 ) have used NH3 and CCS observations to study, among other parameters, 
the gas kinetic temperature in the cores and its variation between isolated and clustered 
environments. 

The above data set of Perseus core observations provides a unique opportunity to carry 
out a statistical, cloud-wide comparison between N2H"'", NH3, and the dust continuum, the 
three most reliable tracers of the dense gas in cores. To this end, we have combined all 
available observations of these tracers in a large set of both pre-stellar and proto-stellar cores, 
and we have carried out a new analysis of the data deriving excitation and column densities 
in an homogeneous manner (checked later against detailed radiative transfer modeling). As 
a result of this work, we present here the first statistically significant comparison between 
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the so-far believed robust tracers of the star-forming core material, and show that while the 
tracers do indeed present similar behavior, significant deviations occur in the densest, most 
evolved core population. 

This paper is organized in the following manner. Section [2] presents the observational 
data used in the analysis. Section |3] compares the physical properties derived from the 
observed N2H+ and NH3 spectra. The chemical properties of the cores are discussed in 
§IU All of the observations are placed in context with core evolution in §0 The major 
conclusions of the paper are summarized in §13 



2. A Coordinated Observational Data Set for Perseus 



As mentioned in the Introduction, the Perseus molecula r cloud has been a mply sur- 
veyed with Spitzer in the mid-infrared (jj0rgensen et al Il2006al : iRebuU et al II2007D and with 



ground-based telescopes at (sub)millimeter wavelength du st emission (IHatchell et al 1 12005 



Enoch et al I l2006t iKirk. Johnstone. Sz Di Francesco II2006I ). At the locations of pre-stellar 
and proto-stellar cores, the Perseus molecular clou d has also been observed in N2H"|" (1-0 ) 
and C^^O (2-1) Jxirk. Johnstone, fc TafallalboO?! ) and NH3 (1,1) JPosolowskv et al Iboosi ). 
As well, near infrared ex tinction mapping using 2MASS sources (for the methodology, see 
Lombardi fc Alves II2OOII) has been used to provide for the large-scale structure of the cloud 
( iKirk. Johnstone, fc Di Francesco 1120061 ). 



2.1. Submillimeter and Mid-Infrared Observations 



Kirk. Johnstone, fc Di Francesco I (|2006[ ) analyzed a 3.5 degree^ region of the Perseus 



molecular cloud using 850 /im data taken with SCUBA at the JCMT. At this wavelength, 
the beam size of the observations is about 15", although smoothing of the map resulted in 
an effective beam of ~20". They identified 58 submillimeter c ores and found that a rnajorit y 
of them could be well fit by stable Bonnor-Ebert spheres ( lEbert I Il955l : iBonnor I Il956l ). 
Comparing the locations of the submillimeter co res with the underlying column density 
in the cloud measured via near infrared extinction, iKirk. Johnstone, fc Di Francesco I (120061 ) 
concluded that the cores are preferentially found in the highest column density regions of the 
molecular cloud. These zones have relatively large r aean densities as well, (n) > 5 x 10^ cm~^ , 
as derived from the extinction maps (see Table 3 in lKirk. Johnstone. Sz Di Francesco II2006I ). 
The majority of the molecular cloud mass, however, was found to exist at low column density 
arguing that only a small fraction of the cloud is participating in the star formation process. 
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A further coi nparison between the 8 50 /im data and Spitzer mid-infrared observations 
was performed by |j0rgensen et al I (120071 ) . In this survey, 72 submilhmeter cores were identi- 
fied (the shght change in number between the two surveys being due to the resolution of the 
reconstructed submilhmeter map and the clump-finding thresholds used to identify objects), 
of which half were identified as harboring protostars. As expected, the submilhmeter cores 
coincident with protostars were found on average to be brighter (more massive) and more 
centrally peaked in appearance. Also, when discernible, the protostars were found to be 
centrally located within the cores. 



A similar map of the Perseus molecular cloud was obtained by lEnoch et al I (120061 ) at 
1.1mm using Bolocam at the Caltech Submilhmeter Observatory. At this wavelength the 
effective beam size of the observations is about 30". The 7.5 degree^ region was inspected 
and 122 cores were identified. In regions of overlap with the 850 fim SCUBA map, the two 
core catalogues are very similar. Due to its larger beam size, the 1.1mm Bolocam data, is 
sensitive to somewhat more extended low-surface brightness sources. 



2.2. N2H+ and CO Observations 



Using the IRAM 30-meter telescope, iKirk. Johnstone, fc Tafalla I (120071 ) simultaneously 
observed C^^O (2-1) and N2H+ (1-0) toward 150 candidate dense cores in the Perseus molec- 
ular cloud. At these transitions, the effective beam sizes are about 11" and 25", respectively. 
For the 89 sources selected by 850 fim emission, 84 % yielded detectable N2H+ emission and 
all were observable in C^^O. The hyperfine structure of the N2H"'" (1-0) emission was utilized 
to fit for the physical properties of the dense gas, including the velocity of the line centroid, 
the line width, the excitation temperature, and the line optical depth. For the C^^O (2-1) 
observations, the line centroid and line width were measured. 



Kirk. Johnstone, fc Tafalla I (120071 ) found that the dense gas associated with the N2H"'" 



pointings displays nearly thermal 



(as determined bv lJ0rgensen et al 



ine w idths, particularly for the subset that appear starless 



20071) . This result is consistent with other surveys of dense 



gas, which used NH3 (IBenson fc Myers Ill989l : iJijina et al Ill999l ). and reinforces the notion 
that cores are supported primarily by thermal pressure, large non-thermal motions having 
disappeared on small scales (< 0.1 pc) and at high densities (> lO'^cm"^). On the other 
hand, the C^^O (2-1) observations revealed that the lower density, un-depleted gas retains 
significant non-thermal motion. Interestingly, the offset between the velocity centroid of the 
C^^O emission and the velocity centroid of the N2H"'' emission was found to be less than the 
sound speed for 90 % of the targets, arguing that the two regions of emission are nevertheless 
coupled (see also § 13. ip . 



- 6 - 



2.3. NH, Observations 



Using tiie GBT, iRosoIowskv et al I ( 120081 ) observed NH3 (1,1) and (2,2) toward 193 
dense core candidates in the Perseus molecular cloud, drawn primarily from the 1.1mm 
and 850 /im continuum surveys. For the observed NH3 transitions, the effective beam size 
is about 30". Ammonia emission was found toward nearly all submillimeter sources and 
the hyperfine structure of the NH3 lines were fit for the observational properties of the 
emitting region, including the velocity of the line centroid, the line width, the excitation 
temperature, and the line optical depths. As well, through comparison of the NH3 (1,1) 
and (2,2) observations, the rotation and kinetic temperature of the gas and the non-thermal 
contribution to the line width was determin ed. For the cores in Pers eus, a typical low kinetic 
temperature of = 11 ±2 K was measured (jRosolowsky et al II2008I ). As found for the N2H"'" 
observations flKirk. Johnstone, fc Tafalla 1120071 ). the NH3 lines are usually quite narrow and 
thermally- dominated . 



2.4. Correlating the Data Sets 

Despite the coordinated approach to data taking inside the Perseus molecular cloud, 
the individual pointing observations in N2H+ and NH3 were not explicitly aimed at the same 
locations on the sky. The N2H"'" observations were taken primarily toward 850 /im SCUBA 
locations, and a subset of 'by eye' extinction locations observed on digitized POSS-II Palomar 
plates. The NH3 observations typically were taken toward peaks in the Bolocam map. In 
general, however, the deviation between the two pointings was significantly less than 25" 
(within the N2H"'"and NH3 beams). 

For the purpose of this study, we collected all observations in N2H"'" and NH3 which 
were less than 25" apart and include them in Table [H Only ten of these sources have offsets 
lar ger than 15" and the mean o ffset is 9". Table [J lists the s ource names from the studies 



by iKirk. Johnstone, fc Tafalla I fl2007l ) and iRosolowsky et al I (|2008[ ). as well as the location 



of the source, and the offset distance between the two surveys. Also presented in Table [T] 
are the integrated line intensities for NH3 (1,1), N2H"'" (1-0), and C^^O (2-1), as well as the 
(sub)millimeter flux from SCUBA (850 /xm) and Bolocam (1.1 mm). The uncertainties in the 
line intensities are typically less then ten percent, while the uncertainty in the (sub) millimeter 
flux is about twenty percent. The (sub) millimeter fluxes are measured at the location of the 
NH3 sources and averaged over 30". 

Tabled] lists 82 individual cores. Of these 82 cores, three have no detection in N2H+ (all 
of these are also found to be weak in NH3). An additional seven cores have poor parameter 
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fits for N2H"'", due to their low optical depth and the tight covariance between excitation 
temperature and optical depth in such a regime (see Appendix One additional core could 
not be fit using the NH3 lines. Thus, of the 82 cores obtained, 71 contain enough information 
to be useful in the detailed comparison of their physical properties. Of t hese, 35 are pre- 



stellar cores and 36 are proto-stellar cores, according to the analysis of |j0rgensen et al 



(I2OO7I ). For these 71 sourc es. Table [2] lists the physical properties derived from fitting to the 
hyperfine components (see lKirk. Johnstone. &: Tafalla 1120071 : iRosolowsky et al 1120081 ). 



3. Physical Properties of the Perseus Cores 

In this section we compare the observationally derived physical parameters for the N2H"'" 
(1-0) and the NH3 (1,1) molecular line transitions. It is worth reminding the reader that 
the effective beam sizes of the two measurements, at their associated telescopes, are 25" 
and 30", respectively. As well, the formal critical densities for thermalizing the emission are 
~ 10^ cm~^ and ~ 10^ cm~^, respect ivel}0. Thus, the ammonia, with its lower critical density 
and somewhat larger beam, should probe to lower density gas if it is present in appreciable 
quantities. 



3.1. Comparison of Centroid Velocities 



Figure [T] plots a histogram of the offset in centroid velocities between the N2H"'" and NH3 
(upper panel) and between the N2H"'" and C^^O (lower pan el). Ten percent of the N2H"'" point- 



ings were found to have multiple velocity components (IKirk. Johnstone, fc Tafalla 



20071 



while considerably fewer of the NH3 observations needed multiple fits (IRosolowsky et al 
20081 ). In this paper, we consider only the closest velocity match between the nitrogen- 
bearing species (and C^'^0) at each position. In this figure (and all subsequent plots), red 
denotes proto-stellar sources while blue denotes pre-stellar sources. Despite the fact that 
the N2H"'" and C^^O observations were taken simultaneously, toward the exact same loca- 
tion on the sky, the N2H"'" and NH3 centroid agreement is significantly stronger. Indeed, 
the mean absolute offset is only 0.07 km s~^, much smaller than the sound speed in the 
gas {cg ~ 0.20 km s~^ for = 11 K) and a bout the accuracy of the No H"*" (1- 0) line rest 



frequency (jPagani et al.l |2009| ) . As noted by iKirk. Johnstone, fc Tafalla I (120071 ) , the mean 



^Note, however, that the critical density is not an exact indicator of the conditions under which emission 
is most efScient. This is especially true for low frequency transitions, such as NH3 (1,1), where stimulated 
emission from the Cosmic Background has a substantial effect on the detailed balance of the energy levels. 
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absolute offset between tlie N2H"'" and C^^O, 0.14kms ^ for pre-stellar cores and 0.17kms ^ 
for proto-stellar cores, is of order tlie sound speed in tlie gas and , interestingly, mucfi smaller 



than the typical C O (2-1) non-thermal line widths (see § l3.2l or lKirk. Johnstone, fc Tafalla 



20071). 



If the N2H"'" and NH3 emission is coming from within the dense core while the C^^O 
emission is produced on larger scales in the material surrounding the dense core, then the 
strong correlation between the centroid velocities of the nitrogen-bearing molecules is ex- 
pected. The critical density for NH3 (1,1) is, however, similar to that required to excite 
the C^^O (2-1) line and thus one might have expected a contamination of the core NH3 
(1,1) measurement by the surrounding C^^O-rich cloud. Unlike Taurus, where the mean 
density in the cloud is low and only the cores reach densities greater than a few xlO^ cm~'^, 
in Perseus the high extinction zones in which the (sub)millimeter cores are found have sig- 
nificant density, (n) > 5 x lO^cm"^, a s derived from the extinction maps (see Table 3 in 



Kirk. Johnstone, fc Di Francesco II2006I ). 



It thus appears that the NH3 (1,1) line is not significantly contaminated by the bulk 
material surrounding the cores in Perseus and that both the observed NH3 (1,1) and N2H"'" 
(1-0) emission is produced within the dense cores. 



3.2. Comparison of Line Widths 



As a second test of the environments in which the various lines are produced. Figure [2] 
shows histograms of the observed line widths (in units of Gaussian a), uncorrected for thermal 
broadeni ng, for N2H"'" (upper panel), NH3 (m iddle panel), and C^^Q ( lower panel). As noted 
by both iKirk. Johnstone, k Tafalla I (12003) and iRosolowskv et al I J2008h . the N2H+ and 
NH3 line widths do not display significant non-thermal motions. This is especially true 
for the measured line widths of the pre-stellar sources (blue histograms), where very few 
measurements fall beyond the sound speed (c^ ~ 0.2kms~^). The C^^O (2-1) line exhibits a 
quite different behavior. In most cases, both pre-stellar and proto-stellar, the measured line 
width is dominated by non-thermal motions. 

As with the line centroid histogram, it is striking how the N2H^ and the NH3 histogram 
measurements agree, in contrast to the comparison with C^^O observations. Despite the 
similar critical densities for NH3 and C^^O, there is no hint of a correlation between these 
two measurements. Again, it would appear that the NH3 (1,1) line is not significantly 
contaminated by the bulk material surrounding the dense cores. 
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3.3. Comparison of Physical Properties 

When determining the column densities and abundances of N2H"'' and NH3 in the next 
section, it is necessary to utilize the observationally fit physical properties to the hyperfine 
structure of each molecule's observed transition. In the preceding sections, the line centroid 
and line width have been investigated. This section compares the key fitting parameters 
source by source. 

In Figure [3l the derived line widths, a^, excitation temperatures, Tex, and line optical 
depths, r, are compared, with the one-to-one ratio shown by the dash-dotted line. Determin- 
ing the optical depths and line widths are fairly straightforward and thus suffer from only 
a moderate uncertainty. The derived excitation temperatures, however, depend strongly on 
the calibration of the instrument and an assumption about the beam-filling nature of the 
emission. Thus the excitation temperatures should be taken as representative values with 
large, at least twenty percent, uncertainties. The top left panel shows a strong correlation 
in the measured line widths of the two nitrogen-bearing species, as discussed earlier in § 13.21 
The top right panel shows that a similar correlation can be found for the derived excitation 
temperatures of the two molecular transitions, although there is a hint of an offset to slightly 
higher, AT ~ 1 K, NH3 temperatures. This result is intriguing given the large difference 
in critical densities of the two species observed. At densities sufficiently large to excite the 
N2H"'' (1-0) transition, nc{^2^^) ~ lO^cm"^, the NH3 (1,1) transition should already be 
thermalized. Unless the gas density is much larger than this critical value, however, the 
N2H"'" should be only sub-thermally excited. That the measured excitation temperatures are 
similar argues that the observed gas in which the N2H"'' and NH3 emission is arising may be 
extremely dense, n ^ ?t,c(N2H"'"). 

The bottom left panel in Figure [3] plots the total line optical depths, integrated over 
all hyperfine components, for the transitions observed. Although there is much scatter, 
the underlying trend is obvious, higher optical depths for NH3 correlate with higher optical 
depths for N2H"'". More importantly, the total optical depths are fairly low - arguing that 
the peak optical depths in any particular hyperfine component are never much larger than 
unity (for NH3 the strongest hyperfine component has a weight of 0.5 while for N2H"'" the 
strongest component has a weight of 0.26). Thus, the integrated line intensities in § 14.11 
should not suffer significantly from optical depth effects. This result has been confirmed 
by comparing the integrated intensity in a single, isolated, hyperfine component against the 
total line intensity. 

Given the expectation that the submillimeter continuum emission at 850 fim is directly 
related to the underlying column density of dust (and by extrapolation the total column 
density of gas - see §|l]and Appendix it is possible to estimate the mean density within 
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each core. In this analysis, we take the 850 /im emission, smoothed to a 30" beam, calculate 
the expected column density of dust and gas, and then divide by twice the measured core 



radius as given in Table 6 of iKirk. Johnstone, fc Tafalla I ( 120071 ) for those cores which have 



radius measures (53 sources). The derived mean density within the core is only approximate 
but may be compared against the derived excitation temperatures to see if there are any 
obvious trends that might be due to sub-thermalization. The bottom right panel in Figure [3] 
shows that for most of the cores the mean density is (ra) ~ 2 x 10^ cm~^. 

The strong correlation in excitation temperature and integrated line intensity, together 
with the equivalence of the kinematic features, suggests that the emission from the two 
nitrogen-bearing molecules is being observed from within a coincident dense region inside 
each core, whether pre-stellar or proto-stellar. In this situation, the true excitation tem- 
peratures are expected to be similar to the derived kinetic temperatures. The fact that 
our measured T^^ value s are often significantly lower than the ~ 1 1 K estimated by 



Rosolowsky et al I (120081 ) is somewhat puzzling, and suggests that our Tex values may have 
been underestimated. This could result from beam dilution effects, if the sources are sig- 
nificantly smaller than the telescope beam, or from additional effects not considered in our 
analysis, like calibration problems or variations of the temperature along the line of sight. 
Detailed modeling of the spatial distribution of Tex and T^ in a selected sample of cores is 
needed to clarify this issue. 



4. Chemical Properties of the Perseus Cores 

In the preceding sections, we have shown that there is a strong correlation between 
the observed properties of the two nitrogen-bearing species, N2H"'" and NH3. The 850 /xm 
submi llimeter flux is also well understood as arising from dust emitting at a temperature T^ ~ 



11 K (IRosolowsky et al II2008I ). In this section we use standard formulae for the conversion 
from observed continuum emission to H2 column density and from observed line parameters 
to N2H"'" and NH3 column densities in order to compare the abundances of these species 
source by source. Although the formulae used here to compute column densities have been 
presented in the literature before, we reproduce them in Appendix |A] for completeness and 
to highlight some misconceptions that often creep into such calculations. 
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4.1. Comparison of Line Strengths 

Given the large number of cores in this sample, 71, it is useful to search for correlations 
in the measured intensities of the molecular lines, as well as against the strength of the 
continuum. Figure H] plots four of these correlations. Where useful, a best-fit linear relation 
is overlaid as a dash-dotted line with the line color denoting the underlying species being fit. 

The top left panel in Figure H] shows the relationship between the 850 /im SCUBA fiux 
(smoothed to a 30" beam) and the C^^O (2-1) emission. There is a possible correlation for 
the pre-stellar cores, although the linear trend is driven by the few C^^O bright outliers. 
For the proto-stellar cores there is no obvious correlation. The top right panel shows the 
relation between the 850 fim SCUBA fiux and the N2H^ (1^0) emission. Again, there is a 
possible correlation for the pre-stellar cores but no obvious trend for the proto-stellar cores. 
The bottom left panel in Figure H] reveals a lack of any correlation between the C^®0 (2-1) 
and the NH3 (1,1) emission, as expected from the discussion in § 13.11 and 13.21 The bottom 
right panel, however, shows a very strong correlation between the NH3 (1,1) intensity and 
the N2H"'" (1-0) intensity, for both pre-stellar and proto-stellar cores. 

For completeness. Figure [5] shows the relationship between the 850 fim SCUBA fiux 
and the 1.1mm Bolocam fiux, both averaged over a 30" beam at the location of the NH3 
measurement. Note that this figure is presented using a logarithmic scale and that the two 
(sub)millimeter measurements correlate exceedingly well. The straight (yellow) line through 
the data is not a fit but rather the expected correlation for dust at = 11 K, and with a 



sub) millimeter emissivity power-law (3 = 2 (for a discussion on /3 see iJohnstone fc Bally 



20061 . and references therein). The typical kinetic temperature of the gas in these sources is 



Tfc ~ IIK (IRosolowsky et al II2008I ). and thus the agreement of the (sub)millimeter measure- 



ments with the 11 K model imply that the gas an d the dust are ne ar thermal equilibrium 



as expected for densities greater than lO'^cm"^ (IGoldsmith I I2OOII ). A few of the bright 
proto-stellar sources have significantly (~ 50 %) higher 850 fim fiuxes than expected which 
may indicate moderate warming toward these locations. In § 14. 21 where we determine column 
densities and abundances, the dust continuum measurements should be reliable to better 
than a factor of 2. 

From the integrated intensity plots in FigurelU it is evident that the C^^O (1-0) emission 
is poorly correlated with the dust emission observed in the (sub)millimeter (especially for 
proto-stellar cores) and uncorrelated with the nitrogen-bearing molecules. For pre-stellar 
cores, there is a possible correlation of the nitrogen-bearing molecules and the dust. 



-VI- 



4.2. Determination of Column Densities 

In general, the (sub)millimeter continuum emission observed from structure within 
molecular clouds is produced by radiating dust grains, and the emission is optically thin. 
Thus, if the dust temperature T^^ and the dust emissivity properties k\ are known, the column 
density of dust is directly computable from the measured (sub)millimeter continuum flux. 
Usually, the dust is assumed to be well coupled to the gas, and the emissivity is given per 
unit gas and dust. The required equation is presented in Appendix lA.il The derived column 
density of H2, A^(H2), toward each pointing is provided in Table [3l assur ning that the dust 



tempe rature is = 11 K, the typical ammonia temperature measured by iRosolowsky et al 



(120081 ). Raising the dust temperature to = 16 K halves the measured column density, 
while a dust temperature of = 26 K is needed to quarter the measured column density 
values. 

Assuming the molecular emission from a given transition is optically thin, the inte- 
grated line strength should be proportional to the column density of the emitting gas (see 
Appendix IA.2p . If the gas is moderately optically thick, and the optical depth can be es- 
timated, correction factors can be utilized to reconstruct the total column density. These 
calculations, however, are extremely dependent on the state of the gas, including the kinetic 
temperature, the density of collision partners, and the equilibrium properties of the molecule. 
For the two nitrogen-bearing species, N2H"'" and NH3, the hyperflne structure of the emission 
provides a reasonable measure of the optical depth, which coupled with the peak intensity 
of the line yields an estimate of the excitation temperature of the transition (in practice a 
more sophisticated flt to the hyperflne structure is utilized). Thus, the number of molecules 
radiating in this transition, along the line of sight, can be deduced. In order to determine the 
total abundance of the molecule, however, the energy partition function must be constructed 
and the level abundances computed. Additionally, for some molecules, including NH3, the 
energy levels are divided between ortho and para forms, which can be treated as independent 
species because they are not connected by normal radiative or collisional transitions. As our 
observations concern only NH3(1,1) and (2,2), in this paper we will only consider the para 
form of ammonia (p-NHs). 

In Table [3l we present the derived column densities for N2H"'' and p-NHs, utilizing 
the formulae presented in Appendix IA.2I The adopted physical properties, cXv, Tgx, and 
r are taken from the physical parameters flt to each spectrum (see Table |2] and § 13. 3p . 
Additionally, determination of the conversion factor, N/Ni, from the column density of the 
observed state, Ni, to total column density, A^, (see Appendix IA.2.2P requires an assumption 
that the energy levels within each molecule are in equipartition at an adopted excitation 
temperature. As discussed in Appendix IA.2.2[ the conversion factor necessary for NH3 (1,1) 
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is relatively independent of the adopted excitation temperature, due to the large energy gap 
between the (1,1) and (2,2) states. We assume the conversion factor for NH3 is 2, which 
is valid to within 10% for temperatures less than 15 K. The conversion factor for N2H"'", 
however, is quite dependent on the assumed excitation temperature used in the partition 
function. In the column density analysis, the derived excitation temperatures for the N2H+ 
(1-0) transition were averaged over all pre-stellar and proto-stellar cores to derive effective 
values of Tex = 5.7 K and T^^ = 6.9 K, respectively. Thus, the conversion to total column 
density for N2H"'' for pre-stellar cores is 2.9, and for proto-stellar cores the conversion is 3.4. 

A test for the validity of the column density ratio measurements between NH.^ and 



N^H"^ was performed using a Monte Carlo radiative transfer co de based on that of iBernes 



(I1979I ) and as discussed in more detail by iTafalla et al I (120041 ). We produced model cores 



with realistic physical conditions for the density distribution, central dust column densities 
similar to those observed for the Perseus core sample, and internal temperatures of 11 K. 
Synthetic spectra were calculated for various input abundances of NH3 and N2H"'", and these 
spectra were reduced using the same procedures as for the Perseus observations. The derived 
column densities and abundances were found to accurately reflect the input values to within 
about 10%. 



In § 13. 3[ the excitation temperatures for both nitrogen-bearing species were seen to be 
similar, suggesting that the density in the emitting region might be larger than the critical 
densities of both molecules. Under this assumption, the true excitation temperature for the 
observed transitions would be ~ Tfc(~ 11 K), as both lines should be thermalized. Deviations 
in the observed Tex from T^ would then be due to beam dilution effects. Under this scenario, 
the conversion to total column density of N2H"'" would be N/Nq = 5.3, or 83% larger than 
the value we have adopted for pre-stellar cores and 56% larger for proto-stellar cores. In 
both cases the uncertainty introduced is less than a factor of 2. The conversion to total 
column density of NH3 would remain 2, however. 



4.3. N2H+ and p-NHa Relative Abundance 

Figure [6] shows the strong correlation between the column density of N2H"'" and P-NH3. 
The tight fit for the two species can be anticipated from the correlation for line intensities 
seen in Figure H] (see also Appendix IA.2p . and the fact that the emission lines are never 
particularly optically thick. The scaling from observed line intensity to column density 
depends both on observing methods (telescope efficiencies and corrections for atmospheric 
attenuation) as well as corrections for the measured physical conditions within the gas (e.g., 
excitation temperatures). Each of these measures has uncertainty associated with it and 
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should dilute any observed underlying column density correlations, as the N2H"'" (1-0) and 
NH3 (1,1) transitions were measured at different times and with different telescopes. As 
well, two different fitting programs were used to determine the physical parameters from 
the hyperfine components. It is more likely that there are correlated uncertainties in the 
integrated intensity of a given species (e.g., an error in the assumed telescope efficiency) and 
thus the scaling for all the column densities of either N2H"'" or p-NHs may be off by about 
20%. 

The top panels of Figure [7] plot the abundance of p-NHs with respect to N2H~'" as a 
function of the total H2 column density (derived from the dust continuum measurements) 
and as a function of the para-ammonia column density. As in all other plots, the blue plus 
signs represent pre-stellar cores and the red diamonds represent proto-stellar cores. The 
scatter in the abundance ratio of these two nitrogen-bearing species is extremely small - less 
than a factor of two about the mean value. As well, there is only a small hint at a difference 
in this ratio for the pre-stellar and proto-stellar cores. For the complete sample of 71 cores 
for which measurements can be made, an average abundance ratio of 22 ± 10 is found for 
p-NHs versus N2H"'". Separating the 35 pre-stellar and 36 proto-stellar cores yields only small 
variations in this abundance ratio, 25 ± 12 and 20 ± 7 respectively (see also Figure E]). 

Given that N2H"'" is a highly reactive m olecular ion, and is kn own to form readily in 



regions where CO depletes (see for example, iBergin fc Tafallal 120071 ). it would appear that 



the formation of p-NHs must follow a similar route to the formation of N2H"'". No discernible 
p-NHs was detected from the moderate density bulk gas surrounding the cores, otherwise, an 
observational correlation with the line widths of the C^*0(2-l) transition should have been 
noted. As well, within the core the line strengths for both nitrogen-bearing species must be 
dominated by dense gas emission in order for the excitation temperatures to be similar (see 
§ 13.31) . Thus, it would appear that both N2H"'" and p-NHs are being produced in the densest 
regions of the cores, in zones where the CO is most likely to be freezing out. 

Perhaps more interestingly, N2H"'" is known to be destroyed by CO and thus in proto- 
stellar cores one might expect that the abundance of N2H"'' should decrease as the internal 
warming evaporates frozen CO back into the gas. No such destruction of p-NHs is postulated, 
however, and thus the lack of any significant change in the relative abundance of these two 
species between pre-stellar and proto-stellar cores suggests that the warming zone is limited 
to a small volume within the core. We return to this discussion on chemistry in §[5l 



- 15 - 



4.4. NHs and H2 Relative Abundance 



The bottom panels of Figure [7| plot the abundance of P-NH3 with respect to H2 as 
a function of the total H2 column density and as a function of the para-ammonia column 
density. For the pre-stellar cores there is only a hint that the abundance ratio is not constant 
at ~ 10^^. The proto-stellar cores, however, show a clear trend of lower ammonia abundance 
in higher column density cores. The trend is reminiscent of the integrated intensity plot 
shown in the top left panel of Figure HI where at high 850 /im flux levels (i.e. high H2 column 
densities) the integrated intensity in the NH3 (1,1) line appears to saturate. As noted in 
§ 13. 3t this is not an optical depth effect; the hyperfine structure of the NH3 does not show 
evidence of reaching the necessarily high optic al depths. A similar effect has also been noted 
for N2H"'" versus H2 for cores in Ophiuchus by iFriesen et al I (120091 ). 



We remind the reader that the H2 column density is derived from dust continuum 
measurements. Consideration of Equation IA2I suggests that significant heating of the proto- 
stellar cores would produce enhanced emission which might be misinterpreted as higher 
column density if the appropriate higher temperature is not used in the analysis. Low 
mass protostars, however, do not have sufficient luminosity to heat large portions of their 
envelope (see for example, |j0rgensen et al Il2006bl : IStamatellos et al 1120071 ). As well, there is 
no evidence for significant warming of the e nvelope in the kinet ic temperature determinations 
derived from the NH3 (1,1) and (2,2) lines (IFoster et al II2009I ). Finally, we note that for the 
brightest 850 /im sources, the calculated under- abundance of p-NHs is greater than an order 
of magnitude, requiring an extreme change in temperature, or other dust properties. We 
return to this issue in §[5l 



Discussion of the Results 



From the observations of N2H+, NH3, C^^O, and dust continuum emission, we have 
uncovered a tight correlation between the abundance of the two nitrogen-bearing species for 
both pre-stellar and proto-stellar cores in the Perseus molecular cloud. No similar correlation 
is found for the nitrogen-bearing species versus CO. As well, only the pre-stellar cores show 
evidence for a constant abundance of these species versus H2. In this section, we consider 
the relevant time scales for chemical and core evolution, in order to better understand the 
possible processes by which the correlations are pr oduced. A more detaile d discussion of 
these processes can be found in the review paper by lBergin fc Tafallal ( 120071 1. 
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5.1. Chemical Evolution 



While the entire chemical pathways to the formation of N2H"'" and NH3 within molecular 
clouds are not fully understood, there is general agreement on the expected steps. For N2H"'", 



(Caselli et al 


1999. 


2002: 


Tafalla et al 


2002) 



has frozen out substantially, the abundance of N2H"'' can increase significantly in the gas 
phase through interactions between N2 and B.^. At the low temperatures expected deep 
within molecular clouds, CO efficiently freezes to dust and thus the relevant time scale for 
CO depletion is the collision time with dust: 



t 



cOjdep 



10^ 



10^ cm' 



yr. 



For NH3 there is no obvious destruction mechanism with CO and thus it is possible 
to have both simultaneously, assuming that there is sufficient time to produce the parent 
product N2 in the bulk cloud. Thus, the absence of NH3 in the extended (C^^O emitting) 
cloud gas suggests that the NH3 formation time scale is longer than the CO formation 
time scale, and comparable to that of core formation and CO freeze out. Indeed, NH3 
(and N2H"'") are considered "late time" molecules due to the long time neede d to activate 



their nitrogen chemistry, which starts with a slow neutral-neutral reaction (e.g. ISuzuki et al 



I992I ). Despite this similar starting point in their production, the different behavior of NH3 
and N2H"'" with respect to the presence of CO in the gas phase makes somewhat surprising 
the strong correlation found between the abundance of these species over the full range of 
core properties seen in Perseus. It should be noted t hat a favored NH3 formation mechanisrn 
starting with the electronic rec ombination o f N2H^ ( IGeppert et al II2004I : lAikawa et al Il2005[ ) 
has been proven inefficient by iMolek et al I ( 120071 ). so a new generation of chemical models 
taking these recent measurements into account is clearly needed. 



5.2. Pre-Stellar Core Evolution 

The formation time for a pre-stellar core depends on the physical processes responsible 
for bringing the material together and is still strongly debated in the community. Before 
the onset of collapse, however, the fastest rate at which the core can be assembled may 
be estimated from the observed physical parameters of the material around the core. The 
maximum infall rate onto the core is 



Mmax ~ 4:11 R^v (n) niriH, 



(2) 
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where mH is the mass of a hydrogen atom, /i is the mean molecular weight, and, for the 
cores in Perseus, (n) ~ 5 x lO'^cm"^ is the mean density in the region surrounding the core, 
i? ~ 5 X 10^^ cm is the typical core size, and v ~ 0.4kms~^ is the typical non-thermal line 
width observed in C^^O (in general one does not expect that this velocity gradient is due 
entirely to infall but using it yields a reasonable upper limit for the infall rate). Thus, even 
if the pre-stellar core is assembled at the maximum rate, a typical Perseus 0.5 Mq core takes 
at least 

W > 1.4 X 10^ ( ) (-—^ X' (^^, ^) yr. (3) 

\5 X 10^^ cm/ v0.4kms~V \5 x lO'^ cm"-^/ 

This time scale is significantly longer than the depletion time for CO within the core, where 
the density is more than an order of magnitude higher. Thus, it would appear that the 
chemical enhancement of N2H"^ should proceed in lock-step with the increase in core mass 
during the pre-stellar core phase as observed (see bottom left panel in Figure [7]). 



5.3. Proto-Stellar Core Evolution 

Once the core becomes sufficiently massive, gravity will dominate over internal thermal 
pressure and the core should collapse on a free-fall time 

,„„.io»(12!£^^l!)"%, (4) 

V ^H2 J 

It is important to note that the time for collapse is longer than the freeze-out time at 
the typical densities inside cores, ~ lO^cm"^, and that the time scale for freeze-out 
drops faster than the collapse time as the core density increases. Thus, even if the core 
were to begin collapse before CO depletion and chemical enrichment of the N2H+ and NH3, 
these nitrogen-bearing species should be substantially enhanced before the collapse becomes 
advanced. 



5.4. The Dichotomy for Dense Gas Tracers in Proto-Stellar Cores 

While the three dense gas tracers, N2H+, NH3, and H2 [traced by the (sub)millimeter 
dust continuum], observed in Perseus appear to mimic one another for pre-stellar cores, 
the trend found between the nitrogen-bearing species and H2 in Figure [7] for proto-stellar 
cores requires a significant drop in the abundance of NH3 and N2H"'" with increasing H2 
column. A similar trend might exist in the pre-stellar core sample but without high column 
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density sources this cannot be confirmed. Had the results been reversed, with the pre- 
stellar cores showing a large range of abundance ratios then one might have appealed to 
a chemical differentiation during core formation. However, assuming that all proto-stellar 
cores began their lives similarly to the observed homogeneous pre-stellar cores, gravitational 
collapse alone (§ 15. 3p does not appear to provide a mechanism for increasing the mass of 
H2 without also increasing the mass of the nitrogen-bearing species. A number of processes 
may be invoked to solve this dichotomy and these are described below. Unfortunately, it 
is not possible to definitively choose which process is most likely, as all have strengths and 
weaknesses. 

One possible solution might invoke the destruction of nitrogen-bearing molecules. In the 
inner region around the proto-stellar core the gas temperature will be raised due to the heat- 
ing from the protostar. Where the gas temperature exceeds ~ 35 — 40 K the CO is observed 



(J0rgensen et al 


2005. 


2002: 


J0rgensen 


2004; 


Dotv et al 


2004) 



and re-enter the gas phase where it will destroy the N2H+. The size of this zone, however, 
is much smaller than the core and thus this effect should be marginal within the single-dish 
observations presented here. For low-mass protostars, the inner warm zone is only ~ 100 AU 
compared with the > 5000 AU core radii, and ~ 7500 AU beam diameters. If the dust evap- 
oration temperature is significantly lower, the evaporation envelope may grow much larger. 
In that case, however, the warmed gas will be rich in CO, providing a strong correlation 
between dust continuum and CO emission for the brightest sources, which is not observed 
in Perseus. Additionally, while the N2H+ abundance will decrease in the inner warm region, 
there is no clear reason why the NH3 abundance should also decrease. 

A second possible mechanism for changing the relative abundances of the nitrogen- 
bearing species versus H2 in proto-stellar cores is for the heavier molecules to freeze-out 
onto the dust particles. Chemical models predict that these species should eventually freeze- 
out, howe ver, there has only been a little observa tional evidence of this effect, especially for 
N2H"'"(see iBergin et al I l2002l : iPagani et al 1120071 . and references therein). In order for this 
scenario to explain the observations in Perseus, the collapsing proto-stellar cores would need 
to be accumulating new material in their outer envelopes, where CO freezes out and N2H"'" 
and NH3 form, while in the interior even these heavy species would be depleting onto dust 
grains. In this manner, the column of H2 would increase over time while the total column 
of the nitrogen-bearing species would remain bounded. 

A third, entirely different explanation for the proto-stellar trend seen in Figure [7] is 
also plausible. While the conversion from emission line strength to total molecular column 
density is relatively straightforward, and apparently validated by the excellent agreement 
between the N2H"'" and NH3 abundances, the conversion from (sub)millimeter dust continuum 
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emission to H2 column density requires knowledge of both the dust temperature and the dust 
emissivity. As mentioned earlier, neither of these values is expected to vary dramatically 
within cores - typical uncertainties in the conversion from observed emission to column 
density for each is a factor of 2, much smaller than the factor of 10 variation seen in the 
observations. We note, however, that it is possible that these uncertainties are correlated and 
that, as the central temperature increases in the core, the dust properties also change. Such 
a scenario could significantly decrease the variation in abundance ratio between the nitrogen- 
bearing species and H2 observed for the Perseus cores. It would, as well, fundamentally affect 
the manner in which proto-stellar envelope masses are determined. 

There is a clear need for the observers and modelers to work together to solve this 
dichotomy between the dense gas tracers if we are to achieve a self-consistent description of 
the proto-stellar environment. 

6. Conclusions 

We have investigated the observed chemical properties of 35 pre-stellar and 36 proto- 
stellar cores in the nearby Perseus molecular cloud. By combining spectroscopic observations 
of N2H"'" (1-0), C^^O (2-1), and NH3 (1,1) along with (sub)millimeter continuum flux mea- 
surements, we are able to determine correlations in the properties of the emitting regions 
for each species, and the abundance ratios between the various chemical tracers. The main 
conclusions from our investigation are: 

(1) The kinematic properties of NH3 and N2H"'" are extremely similar and quite different 
from the kinematic properties of the C^^O molecule, strongly suggesting that the formation 
and destruction of these two nitrogen-bearing species are well coupled. 

(2) For all cores the abundance ratio between the two nitrogen-bearing species is fixed at 
A^(p-NH3)/iV(N2H"'") = 22 ±10. Dividing the cores into pre-stellar and proto-stellar samples 
does not result in significantly different abundance ratios, reinforcing the notion that the 
two nitrogen-bearing species trace the same gas and chemically evolve together. 

(3) For pre-stellar cores the abundance ratio between P-NH3 and H2 is fixed at A^(p- 
NH3)/iV(H2) ~ 10~^, where the H2 column density is derived from the (sub) millimeter 
dust continuum measurements. This reinforces the notion that observations of NH3, N2H"'", 
and (sub)millimeter emission all trace the same dense gas and may be used interchangeably 
when searching for and analyzing pre-stellar cores in molecular clouds. 
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(4) For proto-stellar cores the abundance ratio between the nitrogen-bearing species and 
H2 dechnes from the pre-stellar value as the column of H2 increases, where again the H2 
column density is derived from the (sub)millimeter dust continuum measurements. This 
result suggests that observers should be careful when using a single dense gas tracer (NH3, 
N2H+, or (sub)millimeter emission) to determine the properties of proto-stellar cores. The 
monotonic trend in the observed abundance ratio with H2 column density may indicate a 
simple underlying physical explanation, although decoupling the various possibilities outlined 
in § l5.4l is likely to be tricky and requires further careful observations and detailed modeling. 

7. Acknowledgments 

We thank James Di Francesco, Jon Swift, Paola Caselli, and the anonymous referee for 
helpful discussions. 

Doug Johnstone and Erik Rosolowsky are supported by Natural Sciences and Engineer- 
ing Research Council of Canada (NSERC) Discovery Grants. Helen Kirk is supported by a 
University of Victoria Fellowship. 

We acknowledge the use of data from the following observational facilities: CSO, GET, 
IRAM, JCMT, Palomar, Spitzer , and 2MASS. 



- 21 - 



A. Appendix - Formulae For Determining Column Densities and Abundances 

In this section, we first discuss the underlying relationship between the observed sub- 
millimeter continuum brightness of the cores and the column density of H2. Next we show 
how the integrated line intensities, or physical properties of N2H"'" and NH3 can be used to 
measure column densities. 



A.l. Total Gas Column Density from Submillimeter Continuum Emission 



Assuming a constant dust temperature Td and optically thin conditions, the column 
density of gas, A^H2 1 can be related to the 850 /im submillimeter continuum brightness within 
a beam, 5*850, by 

= S850 [^^bm/i"lH f^SbO Bs5o{Td)]~^ ■ (Al) 

In the above equation, f2bm is the solid angle subtended by the observation, mn is the mass 
of atomic hydrogen, /i = 2.37 is the mean molecular weight, Kgso is the dust opacity per unit 
mass column density (gas plus dust) at 850 /im, and -Bgso is the Planck function evaluated at 
850 /im. Substituting typical values for these quantities, and assuming a 30 arcsecond beam, 
yields 



Nh^ = 3.0 X 10 



22 



850 



1 Jy (30" bm)" 



^850 



0.02 cm-2 g- 



exp 



17K 



- 1 



cm 



(A2) 



Note that the measured column density scales linearly with the submillimeter continuum 
brightness, 5*850, and inversely with the dust opacity, ^850 - Whi le the appropriate value for 
^850 is still uncertain (see for example Ivan der Tak et al Ill999l ). the ranges of values used 
in the literature span only about a factor of 2. In this paper we take ^850 = 0.02 cm^ as a 
fiducial value. 

Only the dust temperature Td enters the equation in a non-linear manner. For the 
sources observed in Perseus, however, the NH3 kinetic temperature measurements suggest 
that the dense gas temperature (which should couple extremely well with the coexistent dust 
temperature) is narrowly scattered around 11 K. 



A. 2. Molecular Species Column Density 



To determine the total column density of a molecule, we first need to calculate the 
column density in the observed transition, A'^; (where the subscript / refers to the lower 
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energy state of the observed transition), and then, through use of an appropriate partition 
function, calculate the total column density of the species, A^. 



A.2.1. Column Density in the Observed Transition 



The column density in the lower energy state / is determined by assuming statistical 
equilibriur n and expressing the o ptical depth in terms of the column density in the lower 
state (e.g., iRohlfs fc Wilson! |2004| ). This relationship is then inverted yielding an expression 
for Ni: 



1 — exp 



1 -1 



(A3) 



where c is the speed of light, h is the Planck constant, k is the Boltzmann constant, Uui is the 
frequency of the {u, I) transition, gi and are the statistical weights, A^i is the Einstein A 
coefficient, Tex is the excitation temperature of the transition, and the frequency- integrated 
line optical depth is J t{u) du. 

For both species in our analysis, p-NHa and the observed transitions have resolved 

hyperfine structure. Hence, r(z^) must reflect the integrated properties of the i hyperflne 
components via 



Tul 



Si exp 



2a,2 



(A4) 



where Tui is the total optical depth in the line as determined by flts to the observed line 
proflle. Si is the weight of the ith hyperflne component deflned such that ^ Sj = 1, z/j is the 
rest frequency of the ith hyperflne transitions and i'lsr is the frequency shift induced by 
the systemic motions of the gas. In Equation IA4t the gas responsible for line formation has 
been assumed to have a Gaussian distribution of line-of-sight motions. The frequency width 
of each hyperflne component is then given by the Doppler relationship: 



0"v 
C 



(A5) 



and is, to high precision, the same for all the hyperflne transitions. Hence, we characterize 
the line with a single rest frequency Uui- Using Equation IA41 the integral in Equation 
can be evaluated yielding 



(A6) 



The relative strengths of the hyperflne components of the transitions allow for a unique 
determination of the optical depth unless the optical depth is very low (see for example. 



- 23 - 



Rosolowsky et al Il2008[ ). Furthermore, the hne widths provide a rehable measure of cXv and 
the hne intensities yield Tex through the observed telescope main beam temperature T^b: 

Tr^y.{y) = Vf [J(Tex) - J(rbg)] [1 - e-^^'^)] , (A7) 

where rjf is the fraction of the telescope beam filled by emission, Tbg = 2.73 K and 

Tui 



J{T) 



exp (Tui/T) 



(At 



Here Tui = {hvui/k), the equivalent temperature of the transition energy. Tmb is related to 
the observed atmosphere corrected antenna temperature T\ by T\ = rj^hi^ui) VfTmb where 
Vmbi^ui) is the main beam efficiency of the telescope at frequency u^i- 



With these observed properties. Equation IA3I can be further simplified to yield 



N, 



4(2 7r)^/^ gi 



1 — exp — 



ul 



T 



(A9) 



It is often convenient to calculate the integrated intensity of the line, as this is a 
robust observational measure. Taking Eqn IA7t converting from frequency to velocity space, 
and integrating over all velocities yields 



lui = [J(Tex) - J(Tbg)] j [I- e-^(^)] dY. 



At low optical depths this converts directly to 

hi = (2 7r)i/2 [J(rex) - J(Tbg)] Tui cr., < 1. 



(AlO) 



(All) 



The last term in the above equation is identical to the last term in Eqn IA9I suggesting that 
we combine the two equations 



lul 



s I 9i J^ul" 



Qu A. 



ul 



[J(Tex) - J(Tbg)] 



1 - exp ( - — 



N, 



(A12) 



where 6{tui, cTv) is a measure of the deviation between the extrapolation of the optically thin 
line intensity equation (Eqn lAlip and the exact calculation. Explicitly, 



5{TuhCr^) 



J [1- e~^(^)] dv 
/r(v)rfv 



(A13) 



This formulation in terms of an escape probability term [6) is useful since it allows the esti- 
mation of column density ratios between two species in terms of their integrated intensities. 
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We append subscripts m and n to indicate quantities for the two species and their respective 
observed transitions: 



'■uLm 



6r, 



1 exp( Tul^m/T^y^^rri) 



uLm 



uLn 



ul,m 



uLn 



(AM) 



L - exp(-T„i,„/Tex,n) 

where Dui^m is a fixed number, dependent only on the physical properties of the molecular 
transition; 

9u,m -^ul,m / \ -I r\ 

Dul,m = 3 • (A15) 

9l,m ^ul,m 

Provided that Tex ^ Tui for both species and Tex,m ~ Tcx,n then the second bracketed 
term in Equation IA14I is approximately unity and the third bracketed term is reasonably 
approximated as T„/,m/^ii«,n, thus rearranging we get 



Ni, 



N, 



Ln 



'■uLm 



(A16) 



Aside from the, possibly large, deviation due to optical depth effects, S{Tui,av), the abun- 
dance ratio of the two states is well represented directly by the integrated intensity. 

An alternate and direct measure of the abundance ratio of these two states can be found 
by taking the ratio of the column densities, via equation IA9t 



N, 



Lm 



N, 



Ln 



exp{-Tul,n/Te: 



1 Gxp( Ty^l,i^/Tg-^„i) 



D 



uLn 



D 



uLm 



Again assuming that Tex ^ T^i for both species and T^.^ ^. 



Tcx,„, this yields 



N, 



N, 



Ln 



Tul,n Dy^i n 



Tul,n C"v,n 



(A17) 



(A18) 



A. 2. 2. Total Column Density Determination 

The observational measures only probe the column density in a single state Ni, which 
can be related to the total column density of the species through the partition function Z. 
Thus, if Zp denotes the number of molecules in state p, 

Z = J2Zp, (A19) 

p 

and the conversion factor from A'^; to is 

N = (^1^ A^,. (A20) 
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For ammonia, there are two terms in the partition function to consider: the distribution 
across the metastable rotational states of the molecule and the distribution between the two 
levels of the inversion transition within each rotational state. Assur ning LTE, the partitio n 
function elements for the metastable rotational states (J = K) are flRohlfs fc Wilson! |2004| ): 



Z, = (2J + l)exp|-A + (A21) 

where B = 298117 MHz and C = 186726 MHz are the rotational constants of the ammonia 



molecule ( Pickett et al 1119981 ) . Care must be taken, however, in calculating the total partition 
function since the ortho- and para- species are not expected to exchange. Thus, there are 
two separate partition functions, one for each of the ortho- and para- states. Since we are 
only treating para-ammonia (p-NHs) in this paper, J 7^ 0, 3, 6, 9... Additionally, the J ^ K 
non- met a- stable states are assumed to carry no weight. 

The full partition function contains a combination of the partition function elements for 
the appropriate metastable states (ortho or para) along with the partition function elements 
of the inversion transition, for which the upper and lower states have equal statistical weight. 
Thus, the full partition function for P-NH3 (J = 1,2,4,5...) is 



^p-NH3 



exp 



(-%^) ] (2J + 1) exp |- A [B J( J + 1) + (C - B)J^] } 



(A22) 

where the first term in square braces accounts for the two inversion levels within each ro- 
tational state. For all inversion transitions, T„;(J) ~ 1 K and therefore the first term in 
the partition function ^ 2 for all relevant core temperatures. The observed p-NHs column 
density, however, corresponds only to the lower inversion level of the J = 1 rotation state. 
When determining 2'i p_nh3 only the first term in the square brace should be included, hence 



Zi. 



P-NH3 



3 exp 



h{B + C) 
kf 



(A23) 



Further, the temperature used in the evaluation of the partition function, T should character- 
ize the level distribution of the ammonia molecules. Since transitions between the metastable 
states are regulated by collisions, this is often equated with the kinetic temperature of the 
gas Tfc or estimated from the observed level populations directly. 



The N2H"'" molecule is a linear rotator, so its partition function is given by 
^N2H+ = J](2J + 1) exp |- A [B J{J + 1)]| . 



(A24) 



where B is the rotational constant: B = 46586.867 MHz ( iPickett et al 
of N2H"'", Zq = 1 corresponding to the J = state. 



19981 ). In the case 



Figure E] plots the conversion factor {N/Ni, where / = 1 for P-NH3 and / = for N2H+) 
as a function of excitation temperature for both P-NH3 and N2H+. For an assumed kinetic 
temperature of = 11 K, the required conversion factors are 2.1 and 5.3, respectively. Care- 
ful consideration of Figure [H] reveals that the conversion for P-NH3 (1,1) is quite insensitive to 
sub-thermal excitation of the meta-stable states (i.e., deviation from LTE). The N2H+ (1-0) 
conversion, however, is quite steep at these excitation temperatures. For example, assuming 
that the effective temperature in the partition function is 6.9 K, similar to the measured exci- 
tation temperatures for proto-stellar cores, yields a conversion for N2H+ of 3.4. An effective 
temperature of 5.7 K, typical for pre-stellar cores in Perseus, yields a conversion of only 2.9. 
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Table 1. Perseus Core Cross- matched Observations 
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Table 1 — Continued 
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0.84 


0.36 


16 


03 


44 


44.2 


+32 


01 


27 


178 


3.8 


48 


1.7 


3.1 


5.9 


1.14 


0.28 


4 


03 


47 


41.8 


+32 


51 


40 


192 


9.3 


49 


6.9 


7.4 


3.9 


0.81 


0.24 
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"From lKlrk. Johnstone, fc Tafalla I 120071 1. 
^From lRosolowsky et al ] I200SI 1. 

'^Separation between iKirk. Johnstone. &c Tafalla I j2007l ) and lRosolowskv et all j 20081 ) pointings. 
■^Within a 30" beam. 



Table 2. Line Properties 



NH3 Properties'' N2H+ Properties'' C^*0 Properties'' 



NH3^ VlSR Cv t T]j Tox Vlsr o"v t Tox Vlsr o-v 

Src (kms-i) (km s-l) (K) (K) (km s'^) (km s"!) (K) (km s"!) (km s"!) 



Pre-Stellar Cores 



8 


4.09 


0, 


.16 


10.9 


11.2 


7.7 


4.05 


0, 


.21 


13.0 


5.7 


4.09 


0.29 


15 


4.65 


0, 


.19 


3.7 


11.2 


6.7 


4.64 


0, 


.20 


4.6 


5.4 


4.55 


0.37 


17 


4.51 


0, 


.15 


14.7 


9.1 


6.8 


4.49 


0, 


.16 


9.1 


5.6 


4.34 


0.37 


29 


5.08 


0, 


.13 


4.7 


10.6 


6.5 


5.09 


0, 


.17 


2.5 


5.7 


4.92 


0.41 


36 


4.71 


0, 


.08 


9.7 


9.1 


6.4 


4.72 


0, 


.08 


3.7 


6.0 


4.89 


0.33 


42 


5.47 


0, 


.14 


3.8 


10.4 


4.9 


5.51 


0, 


.16 


1.0 


9.2 


5.52 


0.20 


47 


8.18 


0, 


.19 


8.2 


11.7 


7.9 


8.19 


0, 


.20 


10.6 


6.0 


8.34 


0.40 


50 


7.21 


0, 


.16 


7.9 


10.5 


6.8 


7.22 


0, 


.16 


6.2 


6.4 


7.28 


0.39 


70 


8.02 


0, 


.36 


3.1 


14.9 


8.7 


7.97 


0, 


.39 


2.4 


11.9 


7.56 


0.80 


72 


8.48 


0, 


.21 


3.3 


12.6 


7.8 


8.48 


0, 


.23 


3.0 


7.6 


8.52 


0.63 


73 


7.76 


0, 


.46 


4.1 


12.3 


8.4 


7.93 


0, 


.32 


5.1 


7.4 


7.75 


0.68 


82.2 


7.51 


0, 


.11 


3.7 


13.4 


5.4 


7.54 


0, 


.07 


10.8 


3.1 


7.53 


0.17 


88.1 


7.41 


0, 


.13 


6.4 


10.3 


6.1 


7.54 


0, 


.21 


7.3 


5.3 


7.53 


0.30 


91 


5.88 


0, 


.14 


7.5 


11.0 


6.9 


5.90 


0, 


.16 


2.5 


8.2 


5.98 


0.28 


95 


6.06 


0, 


.15 


13.8 


10.0 


6.8 


6.08 


0, 


.16 


14.0 


5.1 


6.06 


0.24 


96 


7.81 


0, 


.15 


2.9 


10.5 


5.5 


7.83 


0, 


.16 


8.2 


4.4 


7.78 


0.24 


104 


6.42 


0, 


.15 


7.7 


10.5 


6.1 


6.41 


0, 


.17 


14.6 


4.0 


6.50 


0.42 


105 


6.65 


0, 


.11 


7.6 


9.7 


5.2 


6.65 


0, 


.12 


4.6 


5.8 


6.77 


0.29 


108 


6.80 


0, 


.13 


11.0 


10.2 


6.7 


6.82 


0, 


.16 


9.6 


4.6 


6.57 


0.43 


109 


6.58 


0, 


.23 


1.3 


11.2 


5.2 


6.57 


0, 


.18 


2.7 


4.2 


6.65 


0.32 


111 


6.67 


0, 


.14 


6.2 


10.1 


6.5 


6.65 


0, 


.13 


10.6 


4.6 


6.70 


0.24 


112 


6.58 


0, 


.26 


1.1 


11.5 


5.5 


6.51 


0, 


.21 


9.1 


3.4 


6.63 


0.31 


114 


6.61 


0, 


.19 


9.9 


9.9 


7.1 


6.59 


0, 


.18 


13.5 


4.9 


6.68 


0.34 


127 


6.30 


0, 


.26 


2.1 


11.0 


5.4 


6.22 


0, 


.15 


1.2 


7.1 


6.29 


0.23 


142 


8.47 


0, 


.11 


3.0 


12.2 


5.2 


8.44 


0, 


.10 


4.7 


4.7 


8.39 


0.16 


147 


9.43 


0, 


.10 


9.6 


9.7 


6.6 


9.41 


0, 


.09 


6.7 


5.7 


9.35 


0.14 


148 


9.43 


0, 


.15 


3.6 


17.9 


4.1 


9.40 


0, 


.15 


1.2 


5.1 


9.29 


0.24 


152 


8.82 


0, 


.14 


1.7 


14.0 


4.7 


8.84 


0, 


.09 


2.9 


4.5 


8.74 


0.19 


156 


8.55 


0, 


.35 


1.7 


13.5 


6.1 


8.46 


0, 


.25 


6.1 


4.3 


8.67 


0.42 



00 
to 



Table 2 — Continued 



NH3 Properties'' N2H+ Properties'^ C^*0 Properties'' 





Vlsr 


0"v 


T 










T 


Tex 


Vlsr 


CTv 


Src 


(km s"-') 


(km s~^) 




(K) 


(K) 


(km s~^) 


(km s^-*-) 




(K) 


(km s-i) 


(km s"-') 


157.2 


8.74 


0.14 


3.0 


12.3 


6.6 


8.70 


0.12 


3.0 


5.2 


8.72 


0.24 


169 


8.48 


0.16 


4.2 


11.4 


6.5 


8.55 


0.15 


4.1 


7.1 


8.56 


0.29 


176 


9.97 


0.17 


2.5 


10.8 


5.8 


9.92 


0.22 


4.0 


4.7 


9.59 


0.52 


180 


8.95 


0.12 


4.3 


10.8 


5.8 


8.94 


0.12 


4.9 


5.4 


9.10 


0.36 


189 


10.37 


0.18 


5.1 


10.1 


6.2 


10.37 


0.18 


6.2 


5.3 


10.18 


0.36 


190 


10.13 


0.12 


3.0 


10.1 


5.2 


10.10 


0.13 


2.9 


5.3 


10.06 


0.34 


Proto-StcUar Cores 


7 


4.11 


0.18 


7.0 


12.1 


7.6 


4.08 


0.21 


12.8 


6.1 


4.05 


0.31 


12 


4.53 


0.38 


5.1 


12.6 


8.7 


4.51 


0.35 


6.4 


8.8 


4.35 


0.45 


22 


5.15 


0.15 


6.3 


11.7 


5.7 


5.09 


0.12 


8.9 


5.4 


5.09 


0.19 


31.1 


4.64 


0.17 


3.7 


11.7 


5.8 


4.64 


0.24 


1.9 


8.2 


5.90 


0.19 


32 


4.76 


0.27 


4.4 


13.1 


7.6 


4.72 


0.26 


6.6 


6.3 


4.68 


0.30 


34 


4.97 


0.17 


6.1 


11.4 


8.4 


4.87 


0.20 


6.0 


7.4 


4.91 


0.30 


35 


4.95 


0.21 


3.8 


11.9 


6.8 


5.06 


0.23 


4.6 


6.3 


5.22 


0.45 


40 


7.19 


0.15 


6.9 


10.8 


5.7 


7.21 


0.16 


6.7 


5.6 


7.14 


0.59 


43 


6.82 


0.11 


6.1 


10.6 


5.2 


6.84 


0.13 


3.4 


6.5 


6.86 


0.70 


46 


7.03 


0.15 


7.0 


10.5 


7.4 


7.04 


0.15 


7.1 


6.6 


6.91 


0.66 


48 


7.99 


0.19 


8.8 


11.7 


9.1 


7.98 


0.23 


7.5 


6.8 


8.06 


0.41 


58 


7.48 


0.37 


2.5 


16.5 


8.5 


7.34 


0.24 


4.3 


6.6 


7.62 


0.56 


64 


7.79 


0.31 


2.2 


14.4 


8.7 


7.83 


0.33 


2.2 


9.6 


8.00 


0.65 


65 


7.15 


0.21 


4.5 


12.5 


7.1 


7.16 


0.27 


3.5 


7.6 


7.46 


0.65 


66 


8.01 


0.27 


2.8 


16.4 


10.1 


7.97 


0.25 


4.9 


8.3 


7.73 


0.62 


67 


8.44 


0.30 


3.6 


16.3 


8.8 


8.40 


0.29 


4.7 


10.6 


8.18 


0.64 


68 


7.43 


0.58 


2.3 


16.4 


9.2 


7.74 


0.36 


5.1 


5.7 


7.39 


0.91 


75 


7.32 


0.60 


2.1 


15.1 


7.2 


7.49 


0.41 


5.6 


4.9 


7.30 


0.82 


77 


8.57 


0.22 


3.5 


14.3 


9.3 


8.59 


0.23 


4.3 


8.4 


8.54 


0.60 


78 


7.13 


0.54 


1.9 


14.3 


6.5 


7.15 


0.53 


2.3 


5.2 


6.60 


0.36 


81 


7.49 


0.13 


4.2 


11.8 


6.5 


7.52 


0.14 


4.4 


7.5 


7.55 


0.30 


84 


7.48 


0.15 


3.4 


13.4 


5.1 


7.47 


0.17 


3.8 


5.1 


7.39 


0.44 


87 


7.49 


0.13 


6.6 


10.4 


7.4 


7.50 


0.13 


7.3 


6.8 


7.51 


0.21 



Table 2 — Continued 



NH3 Properties^ N2H+ Properties'' C^**0 Properties'' 





Vlsr 


(Tv 


T 


7k 




Vlsr 


(Tv 


T 




Vlsr 




Src 


(km s~^) 


(km s-l) 




(K) 


(K) 


(km s-i) 


(km s^-') 




(K) 


(km s~^) 


(km s~^) 


89 


8.18 


0.13 


5.5 


10.5 


8.0 


8.15 


0.15 


7.3 


6.3 


8.10 


0.22 


99 


6.99 


0.18 


6.0 


11.4 


7.4 


7.01 


0.25 


2.1 


7.1 


7.08 


0.67 


103 


6.88 


0.24 


5.8 


11.6 


7.7 


6.90 


0.25 


5.3 


7.1 


6.89 


0.41 


118 


6.83 


0.13 


15.9 


9.4 


6.4 


6.84 


0.14 


9.1 


5.3 


6.93 


0.37 


119 


6.41 


0.26 


8.1 


11.5 


7.8 


6.43 


0.24 


13.0 


5.8 


6.65 


0.40 


121 


6.25 


0.33 


5.4 


12.4 


8.3 


6.27 


0.32 


5.9 


7.3 


6.47 


0.51 


123 


6.60 


0.33 


7.2 


11.7 


7.8 


6.66 


0.37 


8.1 


6.5 


6.41 


0.36 


160 


8.64 


0.19 


3.5 


11.7 


8.7 


8.56 


0.18 


4.2 


8.0 


8.52 


0.24 


161 


8.98 


0.22 


4.3 


12.9 


7.7 


9.06 


0.20 


5.1 


6.6 


8.94 


0.34 


162 


8.74 


0.22 


1.7 


14.2 


7.3 


8.74 


0.20 


3.0 


6.3 


8.76 


0.36 


164 


9.00 


0.27 


1.5 


13.0 


7.8 


9.02 


0.28 


2.0 


8.3 


8.82 


0.42 


165 


8.47 


0.18 


5.5 


11.1 


6.3 


8.57 


0.29 


2.2 


6.9 


8.55 


0.26 


192 


10.24 


0.16 


6.0 


11.7 


6.8 


10.24 


0.19 


6.3 


6.5 


10.10 


0.39 



00 



"■From lRosolowskv et "all 1 I2OO8I) . 
''From lKirk. Johnstone, fc Tafallal i20Ql\ ). 
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Table 3. Derived Column Densities for Dense Gas Tracers 



NHs^ 


iV(p-NH3) 


iV(N2H+) 


N{R,) 


Src 


(10^2 cni-2) 


(10^2 cm-2) 


(1020 cm-2) 




Pre- Stellar Cores 


8 


453 


25.3 


363 


15 


167 


8.3 


72 


17 


506 


13.2 


296 


29 


135 


3.9 


81 


36 


181 


2.9 




42 


91 


2.0 


45 


47 


421 


20.8 


314 


50 


293 


10.0 


155 


70 


331 


15.0 


599 


72 


186 


7.7 


287 


73 


547 


18.0 


609 


82.2 


81 


5.3 


296 


88.1 


186 


13.5 


157 


91 


256 


4.9 


106 


95 


478 


19.3 


76 


96 


86 


10.2 


81 


104 


253 


19.0 


141 


105 


157 


5.1 


83 


108 


332 


12.6 


143 


109 


55 


3.9 


91 


111 


197 


11.5 


143 


112 


55 


13.0 


88 


114 


466 


20.3 


197 


127 


106 


1.9 


64 


142 


63 


3.8 


88 


147 


217 


5.9 


76 


148 


81 


1.5 


46 
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Table 3 — Continued 





iV(p-NH3) 


A^(N2H+) 




Src 




(10^2 cm~^) 


(10^0 cm-2) 


152 


40 


2.2 


95 


156 


130 


11.7 


162 


157.2 


93 


3.0 


206 


169 


150 


6.6 


203 


176 


86 


7.2 


160 


180 


106 


5.5 


147 


189 


203 


9.7 


189 


190 


66 


3.3 


23 


Proto-Stellar Cores 


7 


324 


31.3 


613 


12 


573 


33.2 


2100 


22 


186 


11.4 


124 


31.1 


132 


6.4 


151 


32 


316 


20.6 


336 


34 


298 


15.8 


267 


35 


187 


12.8 


233 


40 


210 


11.7 


120 


43 


125 


5.4 


62 


46 


276 


12.4 


179 


48 


524 


21.0 


342 


58 


272 


12.6 


1321 


64 


199 


11.3 


546 


65 


232 


12.8 


114 


66 


260 


17.6 


922 


67 


326 


23.3 


1925 


68 


413 


20.2 


567 


75 


318 


22.9 


3175 


77 


244 


14.2 


812 
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Table 3 — Continued 



Src 


iV(p-NH3) 
(10^2 cm-2) 


A^(N2H+) 
(10^2 cm-2) 


iV(H2) 
(1020 cm-2) 


78 


228 


12.4 


2532 


81 


125 


8.0 


180 




8Q 
oy 


0.41: 


997 


87 


213 


11.6 


158 


89 


189 


12.6 


187 


99 


276 


6.6 


308 


103 


372 


17.1 


718 


118 


464 


12.8 


186 


119 


565 


34.3 


589 


121 


506 


24.8 


968 


123 


650 


35.5 


824 


160 


206 


10.5 


328 


161 


259 


12.5 


896 


162 


90 


7.0 


655 


164 


108 


8.1 


302 


165 


223 


7.9 


253 


192 


236 


14.5 


244 



From 



Rosolowsky et al 



(120081). 
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Fig. 1. — Histograms of the variation in the centroid velocity between the different chemical 
species for the observed cores in Perseus. In both panels the blue lines indicate pre-stellar 
cores while the red lines indicate proto-stellar cores. The top panel compares N2H"'" obser- 
vations with NH3, while the bottom panel compares N2H"'' with C^^O (see text for details). 
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Fig. 2. — Histograms of the measured line widths (in units of Gaussian a) for the observed 
cores in Perseus. In all panels the blue lines indicate pre-stellar cores while the red lines 
indicate proto-stellar cores. The top through bottom panels show the results for N2H"'", 
NH3, and C^^O, respectively. 
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Fig. 3. — Correlations among the physical properties fit to the hyperfine components of 
N2H+ and NH3 for the observed cores in Perseus. In all panels the blue plus signs indicate 
pre-stellar cores while the red diamonds indicate proto-stellar cores. The dash-dotted lines 
denotes a one-to-one correspondence. The top left panel plots the measured line widths 
of N2H"'' against NH3. The top right panel plots the derived excitation temperature, Tgx, 
of N2H"'' against NH3. The bottom left panel plots the derived optical depth, r, of N2H"'" 
against NH3. The bottom right panel plots the derived Tex of N2H"'" against the estimated 
density in the core, ric (see text for details). 
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Fig. 4. — Observed correlations in line intensity, I, and continuum strength, S, for the 
observed cores in Perseus. The symbols have the same meaning as in Figure [31 Where 
useful, a best-fit linear relation is overlaid as a dash-dotted line with the line color denoting 
the underlying species being fit. 
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Fig. 5. — Comparison of the submillimeter continuum flux within a 30 arcsecond beam at 
850 /^m and 1.1mm for the observed cores in Perseus. The predicted correlation, assuming 
a dust temperature = 11 K and a (sub)millimeter emissivity power-law /3 = 2 is shown 
by the yellow line (see text for details) . The symbols have the same meaning as in Figure [3l 
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Fig. 6. — Calculated column densities of P-NH3 versus N2H"'' for the observed cores in 
Perseus. Note that the plot has linear axes. The symbols have the same meaning as in 
Figure [3] and the best linear fits to the data are overlaid as dash-dotted lines with the line 
color denoting the underlying species being fit. The best fit slopes are A^(p-NH3)/A^(N2H"'') 
= 25 ± 12 and 20 ± 7 for the pre-stellar and proto-stellar cores, respectively. 
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Fig. 7. — Comparison of the abundance ratios between N2H"'", P-NH3, and H2 for the observed 
cores in Perseus. The symbols have the same meaning as in Figure [31 The top panels show 
the nitrogen-bearing species abundance ratio against the total column of H2 and P-NH3. 
Here, the central horizontal line marks the mean value of the abundance ratio, 22, while 
the additional horizontal lines indicate a change in the ratio by a factor of 2, up or down. 
The lower panels show the p-NHs to H2 abundance ratio against the total column of H2 and 
P-NH3. Here, the central horizontal line marks the rough value of the abundance ratio, 10~®, 
while the additional horizontal lines indicate a change in the ratio by a factor of 2, up or 
down. Note that the only obvious non-constant ratio is found for the A^(p-NH3)/A^(H2) in 
proto-stellar cores (red diamonds) versus A^(H2) (bottom left panel). 
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Fig. 8. — Conversion factors required to convert the measured column density in an ob- 
served transition to the total column density of the species as a function of temperature (see 
text). The top panel shows the conversion required from A'^i(p-NH3) [the lower state of the 
(1,1) inversion level] to A^(p-NH3). The bottom panel shows the conversion required from 
A/'o(N2H+) [the lower state of the (1-0) transition] to iV(N2H+). 



